The Spectral Decomposition of the Helium atom two-electron configuration in terms 

of Hydrogenic orbitals 
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The two electron configuration in the Helium atom is known to very high precision. Yet, we 
tend to refer to this configuration as a Is f Is 4- singlet, where the designations refer to Hydrogen 
orbitals. The high precision calculations utilize basis sets that are suited for high accuracy and ease 
of calculation, but do not really aid in our understanding of the electron configuration in terms of 
product states of Hydrogen orbitals. Since undergraduate students are generally taught to think of 
Helium, and indeed, the rest of the periodic table, in terms of Hydrogenic orbitals, we present in 
this paper a detailed spectral decomposition of the two electron ground state for Helium in terms of 
these basis states. The Is t Is 4 singlet contributes less than 93% to the ground state configuration, 
with other contributions coming from both bound and continuum Hydrogenic states. 



I. INTRODUCTION 

As early as 1928 Hylleraasi^ recognized that using 
Hydrogenic orbitals to describe the ground state elec- 
tron configuration of the Helium atom was not a good 
idea. In the introduction to his seminal paper he says, 
"It thereby appeared that the use of hydrogen eigenfunc- 
tions, as done by Dr. Biemiiller, leads to erroneous re- 
sults, which obviously has to do with the fact that they 
do not form a complete functional system."— He then goes 
on to introduce a set of coordinates and variational wave 
functions that more accurately describe the two-electron 
ground state of the Helium atomic 

To our knowledge, the results of Dr. Biemiiller were 
never published. Most students nowadays, if asked about 
the electronic ground state configuration of the Helium 
atom, will respond that it is the singlet stated of Is t I s 4: 
as listed, for example, in many periodic tables. Even 
more senior colleagues, while realizing that this simple 
characterization omits interactions, become disturbed at 
the suggestion that other orbitals enter the ground state 
configuration, since this appears, at first glance, to be 
at odds with the 'closed-shell' character of He, and the 
fact that it is an 'inert' element. This misunderstanding 
tends to be reinforced in the way we teach undergradu- 
ate quantum mechanics: since the two-electron problem 
is already very difficult, Helium is often used as an exam- 
ple case study for many approximate methods, including 
perturbation theory and the variational method. Usu- 
ally only corrections to the energy (and not the wave 
function) are described, and even with the variational 
method, often a simple extension of a (singlet) Is t Is 4 
is used. This tends to reinforce the incorrect idea that 
the singlet Is f Is 4 Hydrogenic product wave function 
is the entire pictured 

The solution to this problem attained by Hylleraas was 
to adopt a set of basis functions that more accurately 
capture the electron correlations that exist in the Helium 
two-electron ground state. We want to emphasize here 
that there is no question that this works extremely well. 5 



Nonetheless, physicists and chemists still tend to think of 
atomic electronic configurations in terms of Hydrogenic 
orbitals (the periodic table guides us in this direction), 
and this is how our intuition is formed. Thus, in spite of 
the lesson of Dr. Biemuller's (erroneous) calculation, we 
think it is important to answer the question, 'Just how 
much of the Helium two-electron ground state consists of 
the singlet Is f Is 4 configuration?' and, as an obvious 
follow-up, 'What other Hydrogenic states contribute to 
the ground state configuration?'?. Most (non-expert) col- 
leagues are shocked when the answer requires continuum 
states^ 

Aside from pedagogical value, the requirement of more 
than one electronic configuration has also become a po- 
tentially important ingredient in models of metals and 
superconductors within the field of condensed matter.— ~— 
The electronic ground state of Helium has served as the 
'poster child' for atoms in a solid state environment, 
where the electron occupation can fluctuate. An impor- 
tant ingredient in electron conduction is the fact that 
most lattice models in condensed matter physics focus 
on one set of orbitals from each atomic site; these or- 
bitals overlap to form a conduction band, and the band 
that crosses the Fermi energy, i.e. the energy separating 
occupied from unoccupied bands, is the one band of in- 
terest for low energy excitations. But this reasoning is 
primarily based on single electron ideas, and as a rem- 
edy, this model (tight-binding, or Hiickel) is often gener- 
alized to include electron interactions with one another, 
especially when two electrons occupy the same orbital 
on the same site (the so-called Hubbard model 12 ). The 
Hubbard model, and, in general, Hubbard-like models, 
are the subject of intensive investigation in condensed 
matter*^ However, almost all these models miss the im- 
portant ingredient that the same single electron orbitals 
that are pertinent for single electron occupation are not 
adequate to describe a doubly occupied situation. That 
is, the situation, "when two electrons occupy the same or- 
bital on the same atom" is simply not possible, without a 
modification of those orbitalsJ^ In reality, two electrons 
that find themselves on the same atom will spread out 
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to avoid one another, and therefore inevitably occupy (at 
least partially) other single particle orbitals (like the 2s 
and 3s orbitals in the case of Helium). The present study 
of Helium, while avoiding more complicated ideas like 
Wannier bases and Hartree-Fock orbitals, will emphasize 
this very point. As such, the two-electron configuration 
in the Helium atom serves as a stepping stone between a 
standard quantum mechanical 'textbook' problem and a 
more research-oriented set of problems. 

We first outline the general problem of solving for 
the eigenstates of two-electrons bound to a central nu- 
cleus. This will first be done using 'natural' basis states 
consisting of product states of the bound state hydro- 
genic states. This is presumably the path followed by 
Dr. Biemiiller, but since, to our knowledge, this exer- 
cise has never been published, we provide an outline of 
the method, with most of the details relegated to an ap- 
pendix. As already remarked, this procedure does not 
succeed in converging to the known answers for the He- 
lium atom, for the reason that the continuum states are 
required, and we arrive at the surprising (for some) con- 
clusion that an accurate description of the Helium ground 
state cannot be provided by using just the Hydrogenic 
bound states. 

We thus turn to a more direct approach, where, start- 
ing with the exact wave function, one can compute the 
coefficients of the various basis states of interest. For He- 
lium, while the exact ground state wave function is known 
numerically to many significant digits^ we instead use 
the so-called Hylleraas wave function, 3 consisting of only 
three variational parameters^ in the interest of simplic- 
ity and transparency. This will serve as our 'exact' wave 
function, and, consistent with what we stated above, the 
sum of the magnitudes of the coefficients of the basis 
states consisting of product Hydrogenic orbitals will fall 
short of unity, with the remainder coming from the con- 
tinuum states. 

While much of this discussion is suited for undergrad- 
uates, some of the mathematics is more suited for un- 
dergraduate projects and/or graduate courses. In cither 
case we feel that this approach to the problem provides a 
useful connection between course work and research-type 
problems. We have tried to keep much of the mathemat- 
ics (which we have done in part analytically and in part 
with software packages like Mathematica or Maple) in 
Appendices. 



II. THE HELIUM TWO ELECTRON PROBLEM 

A. Preliminaries 

The energetics of the two electron Helium atom are 
very well known^ Here we concern ourselves with just 
the non-relativistic interactions, so the Hamiltonian gov- 



erning this system is 
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2m 47^0^1 r 2 J 47re r 12 

where r±2 = \fi — ri| is the separation between the two 
electrons, and r\ = \f[\ (r 2 = I^D- Note that we have al- 
ready used an important approximation in writing down 
Eq. ([T]) — we have adopted the Born-Oppenheimer ap- 
proximation, which essentially assumes that the mass of 
the nucleus is infinite. Hence the only degrees of freedom 
are those of the two electrons, each with charge — e and 
mass m. The factor of Z in the middle term is present 
because in general the charge of the nucleus is Ze; obvi- 
ously for Helium Z = 2. It is the last term, representing 
the electron-electron repulsion, that causes the difficulty; 
this is the term which is often ignored (at least concep- 
tually) in a student's first exposure to the Helium atom 
(and the periodic table, for that matter) in his/her un- 
dergraduate education. 

The usual procedure is to first ignore the electron- 
electron repulsion; the problem is then readily solved 
since it now consists of essentially a Hydrogen problem 
that has to accommodate two electrons. The ground 
state solution is then 

i>i{fi,f 2 ) = 0ioo (fl)</>ioo(^), (2) 

which is a product state of two single electron 
solutions; 16 ' 17 

<f> n em(r f ) = Rne(r)Y e m {6,ct>) ) (3) 

where (n£m) are the usual quantum numbers, and R(r) 
[Yfl 71 (9 , cj))] is the radial [angular] part of the wave func- 
tion, which can be written in terms of Associated La- 
guerre polynomials and Spherical Harmonics, respec- 
tively. These functions are tabulated in most texts; for 
reference we write down the ground state result needed 
in Eq. ©: 

Mri) = 2(-f 2 e- Z ^^= (4) 
ao V47T 

where ao is the Bohr radius and it is understood that 
Z = 2 for Helium, but we have left it general in Eq. (0} . 

As it stands, Eq. is deceptively simple; we have left 
out the spin degree of freedom, and therefore we have 
omitted any discussion of parahelium states (symmet- 
ric spatial and antisymmetric spin wave function compo- 
nents) and orthohelium (antisymmetric spatial and sym- 
metric spin wave function components). In the absence 
of spin-orbit coupling these two classes of states are not 
coupled by the Hamiltonian (JXJ) , and one can focus on one 
or the other. Since we are interested in the ground state 
configuration, we will focus on the parahelium states — 
these contain the lowest energy basis state and therefore 
the ground state. This means that all the two particle 
states that we consider below should be understood to 
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include the singlet spin state, 



where 



|x(Si,S'2)) = ^[|t>iU) 2 -U>i|t>2] ) (5) 

where the subscripts 1 and 2 refer to particle 1 and 2J£ 
From this point on spin is no longer included in the dis- 
cussion, but its presence has indeed dictated that we treat 
only symmetric orbital states, of which the state @ is 
one. Other examples are listed in the Appendix, and will 
be used below. 



B. Matrix Mechanics 

Because this approach was so successful in other prob- 
lems, our first line of attack is to utilize matrix mechan- 
ics; we decompose the ground state into a complete set 
of simple well known basis states, 



(6) 



where the a^s are the unknown coefficients. Inserting 
this into the time-independent Schrodinger equation, and 
taking inner products with the same basis states yields 
the eigenvalue equation, 
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47re n - r 2 \ 



(10) 



is the "hard part" , i.e. the electron-electron interaction 
potential. In Eq. © and below we use Z = 2 and 
Eq = h 2 / (2,1710%) ~ 13.606 eV. To consider an actual 
matrix element between any general two electron states, 
we require twelve quantum numbers; hence to avoid pro- 
liferation of indices we use the shorthand i\ = n\t\m\. 
Then 
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((f>ii<Pi2\ H int\<Pi 3 <i>i A ), (11) 



where the last matrix element requires more careful anal- 
ysis; details of the calculation of this matrix element, 
including some examples, are provided in Appendix A. 
There a number of states are listed [see Eqs. (|AlallAlfj) ]. 
and a software package like Maple or Mathematica can be 
tasked with evaluating all the required matrix elements 
up to some cutoff, i max ; the resulting finite matrix can be 
easily diagonalized, and we can obtain the corresponding 
ground state eigenvalue and eigenvector. 



Hij cij = Edi , 
where the matrix elements are given by 



Hn = (^\H\iP 3 



(7) 



(8) 



As already remarked, a complete basis set of Hydrogenic 
orbitals consists of an infinite number of bound states (of 
which Eqs. (lAlaflAlff) are the first few), and an infinite 
continuum of unbound states, for which a much more 
complicated expression is required (though it is essen- 
tially the analytical continuation of the bound states). 
But our tact will be to forge ahead, and simply truncate 
the expansion to include only the low-lying bound states, 
thus excluding a (infinite) number of bound states and all 
the continuum states. Such a truncation scheme worked 
extremely well for the harmonic oscillator^ where, in the 
square well basis used for that problem, only the first ten 
or so states were required to give fully converged results. 

One calculational advantage of this choice of basis 
states is that the one electron parts of the Hamil- 
tonian [all but the last term in Eq. <jTJ)] return 
eigenvalues of the Hydrogen spectrum when operat- 
ing on these states. That is, focusing on a product, 
tpni ,t\ ,mi (rpfini^^jri), which co vers the most general 
case encountered in Eqs. (jAlallAlfj) . we have 

#</W 1 ,m 1 (r!)0n 2 ,£ 2 ,m 2 (^) = 
7 2 



[-Z 2 E Q (— j H j) + -ffint]0ni/i,mi(?'l)0n2,fe,m2( r 2)j 



(9) 



C. Results 

In Fig. 1 we show the ground state energy as a function 
of the inverse of i max , with square symbols, along with 
a guide to the eye; it is clear that the energy has essen- 
tially saturated as i max increases, but to a value much 
higher than the known exact value, which is indicated 
by the large (red) square at the origin. The reason for 
this discrepancy has already been noted: including only 
the bound states is not sufficient — they do not form a 
complete basis set, and the contribution from the contin- 
uum states is significant. At first glance it may seem odd 
that continuum states contribute to the ground state for 
the Helium atom; however, one must keep in mind that 
(i) the right combination of plane waves, for example, 
can indeed describe a very localized state, and (ii) with 
a poor choice of basis states (which, indeed, the Hydro- 
genic states are), the bound states will simply not be 
sufficient to describe the Helium atom. This latter point 
is important, as a truncation in the (infinite) basis for 
the problem in Ref. ((191 ) did not deteriorate the solution 
there, as the set of basis states used in that reference 
were sufficiently 'good' to describe the problem at hand. 
Later we shall demonstrate the role of continuum states 
in this problem. 

Also shown in Fig. 1 is the probability for the singlet 
(Is j" Is I) basis state (right hand scale) as a function of 
1 /imax (blue asterisks) . The horizontal line indicates the 
saturated value, about 91%, indicating that a significant 
fraction (close to 10%) comes from other contributions. 
However, because of the problem mentioned above, this 
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FIG. 1. (color online) Results for the ground state energy 
(filled green squares) as a function of i max (see text and Ap- 
pendix A for a full explanation of i max ). The curve provided 
is merely a guide to the eye. Notice, however, that the ex- 
trapolated result to the origin (as i max — > oo) is well above 
the known exact result, indicated with the large (red) square 
at —79.0 eV. The discrepancy is explained in the text. Also 
shown are blue asterisks, which indicate the probability of 
the (Is t Is 4-) basis state (right hand scale). It saturates 
at approximately (ai) 2 = ct(ioo ioo) ~ 0.91, through which a 
horizontal (blue) line is drawn, indicating that a significant 
amplitude comes from other states. 



value is not reliable, and we will tackle this issue in the 
next section. 

To summarize this section, we have attempted to de- 
scribe the configuration of the two-electron ground state 
in Helium, using only the bound state Hydrogenic or- 
bitals. The reason for doing this is that these orbitals are 
the ones with which we have the most intuition, and they 
are the ones we tend to use to gain a preliminary under- 
standing of the periodic table. Our hope was that, while 
this set of basis states is infinite, a finite set used through 
some truncation scheme would capture the essence of the 
ground state configuration; in reality not only did this not 
work well, it is clear that even if one extrapolates to the 
use of all the Hydrogenic bound states, this incomplete 
basis will not work. 

In fact, the continuum states are required for a proper 
description. The reason for this is that the continuum 
states are able to describe high resolution spatial corre- 
lations that the bound states cannot. The bound states 
become more extended as their quantum numbers in- 
crease, and so any finer scale spatial correlations will have 
to be properly described by the continuum states. This 
conclusion will be reinforced in the next section, where 
we examine projections of the various basis states on a 
very accurate but simple representation of the exact two- 
electron wave function. 



III. PROJECTIONS 
A. Hylleraas wave function: bound states 

One could try to include the continuum states in the 
preceding calculation. To avoid an infinite matrix would 
require a judicious selection of these states, presumably 
based on their energies. In our opinion this procedure is 
fraught with difficulties and various choices are possible, 
so we do not pursue this tact. What we really want is 
a good estimate of the contribution of the 'naive' (Is t 
Is I) basis state to the actual electronic ground state 
of Helium. The contributions from the continuum states 
are of importance only insofar as they influence the value 
a\ ss 0.91 obtained in the previous section. Thus we take 
a different approach, and find that the conclusion of the 
previous section is reinforced quantitatively. 

Using Eq. 0, we can formally rewrite this as 



oo 



m=0 



dp a p \ip c p 



(12) 



where now it is clear that a complete set of states is being 
used; the \ip m ) are those enumerated in Eqs. (jAlallAlfj) 
and beyond, while the |?/>p° nt ) refer to the (as yet un- 
specified) continuum basis states. These states, taken 
together, in fact form a complete, orthonormal set. Mul- 
tiplying on the left by (tp n | and forming the inner product 
therefore projects out the contribution from the n th basis 
state: 



a„ = / i> n ^ 



(13) 



For Eq. (IT31 to be useful we need to know ^cxact! 
as mentioned in the Introduction, this is known numeri- 
cally to many digits precisionj&i£ using a Hylleraas-type 
basis 

* = X) c ij -, / e- fa,/ V ' u l f (14) 

where the Hylleraas coordinates are denned 

s = (r 1 +r 2 )/a , t = {r 2 -r{)/a a , u = r 12 = \fi-f^\/a 

(15) 

and the summation occurs all over non-negative integers 
for j and £, and only even non-negative integers for i. 
This variational basis set has been expanded in a variety 
of ways, as summarized in Ref. (Q), but in what follows 
we take only three terms from Eq. (fT4"|) . with (i,j,£) — 
(0,0,0), (i,j,£) = (2,0,0), and = (0,0,1). Then 

we define the 'Hylleraas wave function' to be the three 
parameter wave function, 



*Hy (Vi\r2) 



2 fkZ 



7r \ ao 



_ — Zks 



[l + 2Z Cl ku + c 2 (2Zkt) 2 ], 

where fc, C\, and c 2 are the three parameters to be deter- 
mined by minimizing the energy, and Z will eventually 
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TABLE I. Results for some overlaps, a; 



i 


basis state 


Oi 


i 1 9 

M 


Tot. Prob. 


1 


100 100 


0.9624 


0.9263 


0.9263 


2 


100 200 


-0.2148 


0.0461 


0.9725 


3 


100 300 


-0.0752 


0.0057 


0.9781 


4 


100 400 


-0.0427 


0.0018 


0.9799 


5 


100 500 


-0.0289 


0.0008 


0.9807 


6 


100 600 


-0.0213 


0.0005 


0.9812 


7 


100 700 


-0.0166 


0.0003 


0.9815 


8 


21-1 211 


0.0260 


0.0007 


0.9822 


9 


210 210 


-0.0184 


0.0003 


0.9825 


10 


200 200 


-0.0146 


0.0002 


0.9827 


11 


200 300 


-0.0090 


0.0001 


0.9828 


12 


100 320 








0.9828 



B. Overlap integrals: continuum states 

For Helium, including the full spectrum of continuum 
two particle states into the problem of computing the 
ground state configuration is very complicated. Noting 
that the most important contributions from the bound 
states arise when one of the electrons is in the one elec- 
tron ground state (see Table 1), we will likewise compute 
only the contributions from the continuum states when 
one of the two electrons is in the one electron (bound) 
ground state. These states, in the singlet configuration, 
are written 

ip P (fi,fi) = -^=(0p(ri)^ioo(r 2 ) + 4>wo{rx)4>p{r2)) (17) 

where only one label (the momentum p = |pj) is required 
because the other electron is in the Is state so the con- 
tinuum state also has I = m = 0, and we have written 
the right-hand-side of Eq. (fTT)) as depending only on the 
radial coordinates, r\ and r 2 . The single particle state 
4> p {r) has a radial part R P (r) given by 



be set to equal 2. Eq. (|16|) will serve as our 'exact' wave 
function. In reality, minimization of the energy with ^Hy 
yields E^y ~ —78.979 eV, with variational parameters 
ci w 0.0803 and c 2 0.0099, and k « 0.908. Evaluating 
the necessary integrals and obtaining these results is well 
documentedf^i and so these are simply summarized in 
Appendix B. The attained energy is within 0.05% of the 
exact result (—79.014 eV); we consider this sufficiently 
close to justify our adoption for the present purposes of 
-^=| v E , Hy) as the 'exact' normalized Helium wave func- 
tion, with N given below. The calculation of the overlap 
integrals for the bound states can be done for the general 
case — some details are provided in Appendix B, along 
with a simple example. Our results are summarized in 
Table 1. 



It is clear that the largest contributions arise from 
states in which one of the electrons is in the 4>wo (r) state, 
i.e. the single electron ground state. In particular, more 
than 92% comes from the Is t Is I singlet, but a sizeable 
contribution comes from states other than this one. Fur- 
ther (small) contributions arise from states not listed. 
While the total probability (5th column in Table 1) is 
close to unity, inclusion of the remaining bound states 
(not shown) still gives a total probability that falls short 
of 0.99. The remaining probability arises from continuum 
states, to which we turn in the next subsection. 



R p (r) 



27r£fi 



1 — e pa o 



-ipr 



M(l 



i ,2,2ipr), 

pao 



where we have written this for general Z (here we need 
only Z = 2), and 



M(a,b,z) 



E(a) m z m 
m—0 v ' 



(19) 



is the so-called Kummer function^ and (a) rn = a(a + 
l)(a + 2)...(a + m— 1) is the so-called Pochhammer sym- 
bol. Note that (a) = 1- Equation (fT5|) is essentially the 
analytical continuation of the radial bound state wave 
functions. The standard^ normalization condition for 
the continuum states, 



f 



drr 2 R p ,(r)R p {r)=8(p l -p), 



(20) 



determines the coefficient in Eq. ([18)) . Just as for the 
bound states we require 



»p- j d 6 n J d 6 r 2 tl) p {r 1 ,r2)-j=^Hy(fi,f 2 ). (21) 

Because of symmetry the contributions from the two 
terms in Eq. (fl7|) are identical; the integrals involved 
in Eq. (|21[) can be done, either with Mathematica or 
by hand. Some detail is provided in Appendix C. With 
y = pao I Z the result is 
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ci [I6fc/i + 4fc(fc + l) 2 / 3 - Ak(k + 1) J 2 - 16fcJi] + 

c 2 [96fc 2 / 2 + 8k 2 (k + ifh - 48k 2 (k + 1)J 3 ] 1 (22) 



where the expression for 7„ and J n are provided in Ap- 
pendix C. Then the total contribution for these contin- 
uum states, 

2 f°° 

fcont = - / dp\a p \ 2 , (23) 

is evaluated numerically; the result is a further contribu- 
tion of Pcont ~ 0.0117, which brings the total probabil- 
ity from the bound states in Table 1 and the continuum 
states shown in Eq. (|17|) to ~ 0.995. The remaining prob- 
ability required to reach unity comes from bound state 
contributions not included in Table 1 and from contin- 
uum states beyond those not considered in Eq. ([T7]) . 



IV. SUMMARY 

The primary purpose of this paper was to demonstrate 
the degree to which the electron configuration in the He- 
lium atom is not simply the singlet Is J Is J, state. In 
fact we have shown, through two methods, diagonaliza- 
tion in a particular basis set, and through projection on 
this same basis, that approximately 8% of the wave func- 
tion is not the Is t Is I singlet. As already noted, this 
has a very large effect on the energy. While we mentioned 
in the introduction that this fact may be of interest in 
current research in strongly correlated electron systems, 
including superconductors, the emphasis here has been 
on pedagogical aspects. Further exploration of the con- 
sequences of this electron configuration in a solid can be 
the topic of specialized student research projects. 

The diagonalization of the problem in a truncated basis 
was, in fact, not successful at reproducing the ground 
state energy; in being unsuccessful, this calculation has 
served to highlight important pedagogical points, and it 
is partly for this reason that we have included a detailed 
analysis of the problem here. The two-electron problem 
in Helium served to highlight that the bound Hydrogenic 
states do not form a complete basis set, and while one 
might have thought that the low-lying states would be 
sufficient to accurately describe the ground state, this 
work shows that this supposition is incorrect. 

Setting up the electron configuration in Helium as a 
matrix diagonalization problem also serves to provide a 
concrete example of matrix mechanics, which, for un- 
dergraduates at least, is often introduced only in a for- 



I 

mal sense, with an abstract-only exposition of Hilbert 
space and expansion in basis functions, etc. This prob- 
lem also serves to provide a first exposure to a realis- 
tic 'many-body' problem, and how one would construct 
many particle wave functions (here, two electrons is the 
first stepping stone in this direction beyond the single 
particle problem to which undergraduates are normally 
exposed) . 

The second method we presented, projection of an ac- 
curate wave function onto an orthonormal basis set of 
'known' and well-understood wave functions, requires a 
little more sophistication, because Hylleraas wave func- 
tions are usually only introduced in graduate studies. 
Nonetheless, the three-parameter Hylleraas wave func- 
tion given in Eq. (fT6|) is sufficiently simple to be suitable 
for undergraduates. The ensuing algebra for the bound 
states is also readily accessible, particularly if projec- 
tions are computed on a case-by-case basis. For example, 
evaluation of Eq. (|B11[) requires knowledge only of ele- 
mentary functions and standard integrations. We have 
included projections involving the continuum states as 
well, and these require more advanced mathematics, and 
certainly can be skipped at the undergraduate level. 

The topic of even two electron correlations has been 
traditionally confined to more advanced studies at the 
graduate level. Part of the reason for this is that there are 
better and more accurate methods available after a pre- 
liminary introduction to many-body methods has been 
assimilated by the student. However, most students see 
the electronic structure of atoms in the periodic table in 
general, and the Helium atom in particular, in terms of 
Hydrogen orbitals, and therefore it is desirable that a de- 
scription be provided in terms of these same orbitals, as 
we have presented in this paper. 
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Appendix A: Matrix mechanics for Helium example: 
1. General considerations 



The first 6 two particle states that need to be consid- 
ered [dropping the spin part given in Eq. ([5])] are, for 



V2 



0ioo(n)0ioo(f*2 
[0ioo(?"~i)02oo(r2) + ^loo^D^ooW) 



V2 



[(f>iaa(ri)(f> 2 i-i(f 2 ) + 0ioo(n2)02i-i(rl) 



1 
1 

71 



[01Oo(n)021o(r2) + 0ioo(?"2)02io(ri) 

[01Oo(?'l)021l(^) + 0ioo(?*2)021l(n) 
0200 (?"!) </>200 ipl 



m,\ + m2 = 
mi + m 2 = 



uii + m2 



-1 



mi + ?ri2 = 

mi + ?7i,2 = 1 
mi + m,2 = 



(Ala) 
(Alb) 

(Ale) 

(Aid) 

(Ale) 
(Alf) 



Note that all these states are symmetric, and that we 
have indicated the quantum number mi + m 2 for each 
two-particle state. In any state in which this is not zero, 
that state will not contribute to the ground state. This 
selection rule will be derived below, along with some 
other rules that eliminate more of these states. These 
rules provide a big simplification, and to the degree that 
we anticipate that maybe the [nim) = (600) doesn't 
contribute very much to the ground state, the number 



of states that may be required might be anticipated to 
be small. Unfortunately, this is not the the 
continuum states contribute as well; this fact ultimately 
thwarts this attempt to attain a spectral decomposition 
in terms of these basis states^ Nonetheless, we proceed 
along this line of investigation, first to demonstrate that 
this is the case, and secondly to establish these helpful 
selection rules. 

To evaluate the last matrix element in Eq. (fTTj) we 
introduce the expansio n 23 ' 24 



n - r 2 \ 



E 



(k - |m fe |) 



■ r < p\ m k\ 



{cos9 1 )Pl mkl {cos9 2 )e lmk( - 



(A2) 



where r< = r\ and r> = r 2 if T\ < r 2 and vice versa where e m = (— l) m for m > 0, and e m = 1 for m < 0. If 
for n > r 2 . The P™ fc are the associated Legendre poly- we substitute into the last term of Eq. (fTTj) ) we obtain 
normals, and the angles correspond to the spherical co- 
ordinates for each of the vectors, f{ and f 2 . Eq. (|A2j) 
should become familiar to students through problems in 
Electromagnetism as well as in Quantum Mechanics. We 
can expand each of the individual Hydrogenic states (see 
Eq. ©) and use 



V 



' ( 2£+ l)( £ _ H)! 

4tt (e+\m\)\ e Fi { - C0S&h 

(A3) 
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^■> i2 \H int \4>iM = £ £ §-4^(~i) (mi+|mil+m2+l " l2l+m3+|m3l+m4+|m4l)/2 



' (2Zi + lWi-KI)! (2l 2 + l)(l 2 -\m 2 \)\ / (2l 3 + l)(l 3 - \m 3 \)\ / (2Z 4 + 1)(Z 4 - |m 4 |)! 
(Zi + |mi|)! V (/ 2 + |m 2 |)! V (f 3 + |m 3 |)! V (k + |m 4 |)! 



/ Rn 1 .,h( r i) R n2,h( r 2)Rn 3 .j 3 (r 1 )R niM (r 2 ) n ^ T rlrldr 1 dr2 
lo Jo r% +L 

i r 



X 2 



/ P; 1 ™ 1 \cos6 1 ) F ; ' 3 '™ 3 1 (cos^i )P| mfc 1 (cos^i )sind 1 d8 1 
Jo 1 



x r P 'r 2 ' ( cos Wf 41 (co^ 2 )pl mfc| {cos6 2 )sine 2 d6 2 
2 Jo 



fc 

27T J 2ir J 



X 



2. Simplifications and Selection Rules 

The last two integrals in Eq. (|A4[) require to be 
fixed; moreover, compatibility between the two leads to a 
condition on quantum numbers of the states that lead to 
a non-zero expectation value of the interaction potential: 

mi + m 2 = m 3 + m 4 . (A5) 

Inspection of the states (jAlal - lAlfjl and those beyond 
shows that many states (e.g. lAlcl and lAle[) do not con- 
tribute to the ground state, and hence can be discarded 
from further discussion. 

Furthermore, Gaunt was able to evaluate the 9\ and 2 
integrals analytically. 2 ^ We have actually found it simpler 



to evaluate the integrals as they are (either numerically 
or analytically with the aid of Maple or Mathematica) ; 
nonetheless Gaunt 's formula leads to 2 ^ the so-called tri- 
angular condition for the angular momenta. In the first 
(6*i ) integration in Eq. (|A4[) . for example, this requires 

£ 1 + 4>fc>|£ 1 -£ 3 |. (A6) 

Thus, the &-sum in Eq. (|A4[) is truncated at k maXl where 

'Mn ax 

= min{l\ + £ 3 , l 2 + £4). (A7) 

It is conventiona l 24 ' 26 to introduce coefficients defined as 
follows: 



c k {£rrf£'m') = (_i)( m +M+ m '+l m 'l+ m - m '+l m - m 'l)/ 2 



(k-\m-m'\)\ (21 + l)(l- H)! {2C + !){£'- \m'\)\ 
y k+\m-m'\)\ y (£+\m\)l \] (£' + \m'\)l 

x~ f_' P\ m \»)P^\»)P\r- m X^. (A8) 

In terms of these coefficients the interaction matrix element can be written as 

e 2 X 

(fa^i^Hint^i^is) = — i mi+ra2:m3+m4 } y c k (£imi; £ 3 m 3 )c k (£ i m i ; £ 2 m 2 )R k (n 1 £ 1 ;n 2 £ 2 ; n 3 £ 3 ; n 4 l 4 ), (A9) 
where 



k=0 



poo poo fc 

R k (n 1 £ 1 ;n 2 £ 2 ;n 3 e 3 ;n 4 £ 4 ) = J J P„ 1 / 1 (ri)P„ 2!fe (r 2 )i?„ 3!fe (r 1 )i?„ 4i£4 (r 2 )-^ T r 2 r 2 dr 1 dr 2 (A10) 



This last double integration can be readily done by hand Mathematica or Maple, since the radial wave functions 
(though it is tedious) or can be done (symbolically) with can be readily expressed in terms of associated Laguerre 
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polynomials. 17 Also note the reversed order of the argu- 
ments in the second c k in Eq. (| A9|) ; this is important 
since 



c k (£m;£'m') = (-l)(™'- m ) c * (£'m' ';£m). 



(All) 



Inspection of the states (jAlallAlfj) and those beyond 
not already omitted by the condition (|A5[) indicates that 
three distinct possibilities remain (e.g. when (niiiirii) = 
(riiiim-i) and (n-ii-im^) — (n^^m^), etc.) but all of these 
can be handled through Eq. (|A9I) . 

One other selection rule is present in these results, 
though not readily apparent. In Eq. (|A9[) the same value 
of the index k must work for both c k coefficients. These 
coefficients have been tabulated, for example, in the texts 
by Slater^ where it is clear that two even or two odd 
angular momenta (referring to £ and £') couple to one 
another only through even values of k, while an even and 
an odd angular momentum couple to one another only 
through an odd value of k. For example I = and £' = 
(ss) results in a non-zero c k coefficient only if k = 0, 
while I = and £' = 2 (sd) yields a non-zero c k coeffi- 
cient only if k = 2; the case £ = 1 and £' = 1 (pp) has a 
non-zero c k coefficient if k = or k = 2. In contrast the 
sp, sf, and pd coefficients are non-zero for k = 1, k = 3, 
and k = 1,3, respectively. This means that £\ +£3 has to 
have the same parity as £2 + £4- Since the parity of the 
Is T Is 1 state is even, then, for example, we can further 
discard state (|Ald[) from the list, leaving only 3 of the 
original 6 states listed. 



3. An example 

By way of example we quote some steps for the first 
state (not listed) that does not utilize single particle 
s states, 

1pu(fl,r2) = [021-l(f"l)021l(7'2)+'/'21-l(r2)(/'21l(? , l)] . 

To compute the matrix element that couples the first 
(and primary) basis state with this one, we need to eval- 
uate the overlap integral, 



#1,14 = (V'll-HintlV'14)- 



(A13) 



This consists of two overlaps (Eq. (IA12[) has two terms) , 
but, through a change of variables f% O fi, these are 



readily seen to equal one another. We are thus left with 

#1,14 = V2(</>ioo01Oo|#int|</>21-1021l)- (A14) 



Eq. (TA6j) tells us that I > k > 1, i.e. only k = I 
need be considered. Then a straightforward evaluation 
of c^OO; 1 - 1) and c x (ll;00) (as required in Eq. (fS9|) ) 
gives us -1/Vo for the first and +l/\/3 for the second 
(using Eq. (|A11|I 'I. Next only one double radial integral 
(Eq. (|A10[) for k = 1) is required, and a straightforward 
evaluation, readily done by hand, gives 



^(10; 10; 21; 21) = 



112 Z 
2187 a^' 



(A15) 



Combining this with the c 1 's and the V2 from Eq. (|A14p 
gives 



Hi, 



/tt 448 „ 
-^6561*' 



(A16) 



where E, 







ti 1 



13.606 eV is adopted as our unit of 

energy, and we have used a Q = ^fl 2 -^- 

In practice we have written a program in Maple to per- 
form these tasks, for the bound basis states, up to some 
cutoff. That is, all the matrix elements, Hij, up to some 
cutoff, imax, are evaluated, and then this matrix is diag- 
onalized. The cutoff i ma x is defined as the 'n ! quantum 
number up to which all states are included. For example, 
if *max = 2, then only 5 basis states are considered. The 
rest of the states either do not contribute to the ground 
state, or, if one of the Hydrogenic single particle states 
has n = 3, only contributes to the next shell (i max = 3) 
and beyond. Matrix elements for the case « m ax = 2 are 
tabulated in the next section. 



4. Matrix equation for / max = 2 

If we restrict basis states in Section II to those with n = 
2 or less, only 5 two particle states need to be considered; 
referring to Eqs. (|AlallAlf| . these are tp2, tpe, ^14, 
and i/?i6 , where ipi^ is given in Eq. (TAT2|) . and 



4>w(ri,fi) = 0210(^1)^210(^2). 



(A17) 



Writing the wave function in terms of these wave func- 
tions alone gives rise to a 5 x 5 matrix diagonalization 
problem (see Eq. ([?]))• The resulting equation is: 



_L_L 

2 

32768 
64827 
64 
729 
448 
' 656 
448 
6561 



\V2 



32768 
64827 
2969 
729 
4096 
84375 
2048 
28125 V 
2048 



2X125 



64 
729 
4096 
84375 
179 
128 

15 /o 
128 VZ 
15 
128 



2048 ^2 



28125 

128 VZ 
47 
40 

— 22- a/2 

640 v z 



448 
6561 
2048 
28125 
15 
128 
27 p 
" 640 v< 
779 
640 





ai 




ai 




a 2 


E 


a-2 




a 7 




a 7 




ai 4 




«14 




. ai 6 _ 




. «16 . 



(A18) 
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where Eq « 13.606 and the subscripts on the coefficients 
correspond to the labels in the wave functions. The re- 
sulting ground state energy is E\ = —77.13 eV, and the 
eigenvector has components 



°100,100 
<2l00,200 
^200,200 
d21-l,211 
1210,210 



0.9520 
-0.3040 
-0.0146 

0.0267 
-0.0188. 



(A19) 



basis state dominates, and small ad- 



Clearlythe ^(ioo.ioo) 
justments occur as the number of basis states increases. 
Nonetheless, almost 10% of the wave function is com- 
prised of components beyond the Is f Is I state. The 
matrix construction and diagonalization indicated here is 
repeated, with the aid of the software package MAPLE, 
for increasing values of I max , and the results from these 
calculations are reported in the text. 



Appendix B: Hylleraas wave function 
1. The expectation value of the energy 



For reference, the expectation value of the energy, us- 
ing the normalized Hylleraas wave function — 3=|\]/jj y ) 
with I^Hy) given by Eq. (fT6|) . is 



Eiiy — — 



AM 



2mal N 



{kZf 



(Bl) 



where k can be determined analytically in terms of c\ 
and C2 [see Eqs. (|B4[) . below], 



kZ 



ZL-L' 
2M ' 



and the unit of energy is the Rydberg, 



1 Ryd sa 13.606 eV. 



2maQ 



(B2) 



(B3) 



The parameters L, L' , M, and N, corresponding to dif- 
ferent parts of the Hamiltonian, are given by 

L =4 + 30ci + 48c 2 + 72c 2 + 280cic 2 + 576^ 
L' =5/4 + 8ci + 9c 2 + (35/2)c 2 + 48cic 2 + 78c 2 



M =2 + (25/2)ci + 24c 2 + 32c? + 146cic 2 + 480c 2 . 



N =4 + 35ci 



48c 2 + 96c 2 + 308cic 2 + 576c 2 .. 



(B4) 



Aside from units (most treatments use so-called chem- 
istry units) these all agree with results in the 
literatures^ Use of the optimally determined param- 
eters, ci w 0.0803 and c 2 w 0.0099, and k « 0.908 deter- 
mines the energy via Eq. (IB1|) . 



2. Overlap integrals: bound states 

With a very accurate wave function in hand, we can 
simply utilize Eq. (fl~3|) to determine the probability of 
each basis state in the ground state wave function. The 
general Hydrogenic bound state can be written as 



where 



Gn,l — 




£)\} 3 \na J 



(B5) 



(B6) 



and the required overlap integrals consist only of terms 
of the form 



1 



N 



Hy 



(B7) 



In the u,s,t coordinates defined above, and using the 
volume element, 



dr = — (s 2 - t 2 )u dsdtdu shi0id0i 
8 



(B8) 



many of these integrals are straightforward to calculate 
analytically — an example will be shown below. How- 
ever, it is of interest to compute these overlaps for gen- 
eral quantum numbers, in particular to check on selec- 
tion rules, and to examine the convergent behaviour for 
large m and n 2 . To this end we use the same expansion, 
Eq. (|A2[) used earlier to evaluate matrix elements of the 
Hamiltonian, and then evaluate Eq. (|B7[) . The details of 
the derivation for a general matrix element are a little too 
cumbersome to include here; all of the overlaps involv- 
ing bound state wave functions can be expressed in terms 
of Hypergeometric functions, which are readily evaluated 
using Mathematica, and with these we can readily sum 
the contributions from the bound states. 

For students, however, it is more instructive to com- 
pute this overlap 'by hand' for a few of the most relevant 
states. We outline the procedure for the most important 
state, ii\{f{,fi) = 0ioo (rl) (j> 100 (rjj). Then 



ai 



d 3 n I d i r 2 c>ioo (rl) 0ioo (»*2)-^=*Hy(n,rl), 



N 



. {B9 ) 

and the only angular dependence occurs in the |rj — r 2 | 
term in the Hylleraas wave function. This requires knowl- 
edge of the angle between the two vectors, rj and r 2 , 
which we can call 0i 2 ; this is given as^ 

cos0i 2 = cos 0i cos0 2 + sin0i sin0 2 cos (0 2 — <f>\), (B10) 

where fl = (n, 9i, 0,-), for i = 1, 2 are the spherical coor- 
dinates for each vector. A simpler trick is to line up the 
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z 2 axis with f±; then #12 = 6 2 , and all six integrations 
are then straightforward. The result is 

32 fc 3 r k , k ,2! 

a x = -=- -<U + 35 Cl - + 96c 2 (- ) }, 

V^v(fc + i) 6 l fc + i v fc + i ; J 

w 0.9624, (Bll) 

and therefore a\ sa 0.9262. The numerical values follow 
upon substitution of the optimal values of fc, c\, and c 2 . 
Similar calculations can be performed for the other over- 
lap coefficients; some results are tabulated in Table 1, 
where the constituents of the basis states, suitably sym- 
metrized, are listed. Note that the label provided under 
the 'basis state' column identifies the two single particle 



wave functions involved, and the basis state is a singlet 
state, and therefore the spatial part is symmetrized, as 
enumerated in Eqs. (|AlallAlfj) . 

Appendix C: Evaluation of overlap integrals 

We provide some of the details for the evaluation of 
Eq. (l2~Tj) . Because of the symmetry in the singlet state , 
and using the definition 




we require the integral 



f fc fc 

A p = I d 3 n / d 3 r 2 p (r 1 )</) 1 oo(r2)e"^ (ri+r2) {l + 2Zc 1 — |rl-r2|+c 2 [2Z— (n-r 2 )] 2 }, (C2) 
J ao ao 



where 



1 y / 27T— 7 

6p(n) = -= — \- , e- ipr M(l + i ,2,2ipr), (C3) 



and M(a,b,z) is the Kummer function. Note that 



/ d6 1 sm6 1 de 2 sme 2 dfa d<t> 2 [1] = (4tt) 2 , (C4) 
Jo Jo Jo Jo 



10 Jo 

while 

c27T r 2ir 1 „2 



/ d0 1 sm0 1 d9 2 sm0 2 dfa dnh\f{ - r 2 \ = (47r) 2 (r> + --^), (C5) 
Jo Jo Jo Jo r > 



where r< (r> ) refers the smaller (larger) of n or r 2 . The 
first of these is trivial, while the second follows most read- 
ily by using the trick mentioned following Eq. (|B10|) . 

Two integrations remain, and it is efficient to switch to 
dimensionless variables, x\ = Zri/cto, x 2 = Zr 2 /ao, and 
y = pao/Z. Then, since the Kummer function depends 
only on x\, the x 2 integration is elementary, and leaves 
behind polynomials in x\. One arrives at the expression 
given by Eq. (f2"2"j) . with the definitions 



I n (y) = B n (y, k + iy) and 
J n (y) = B n {y,2k + l + iy) 7 



(C6) 



where 



Bjy,z)= dx x n e- zx M(l + -,2,2iyx). (C7) 
Jo V 

Now an expansion of the Kummer function (following Eq. 
(HHJ)) allows us to do the integral and obtain an infinite 
summation which can be recognized as a Hypergeometric 



function, 21 



B n (y,z) = -^F[l + -,n + l;2;^], (C8) 
z™ +i y z 



with n — 1, 2, 3, 4 needed. The standard definition of the 
Hypergeometric function uses the Gauss hypergeometric 
series with circle of convergence in the unit circle \z\ = 
1 



F(a, b; c; z) 



(a) m (b)m Z r 



00 



ml 



(C9) 



where the (a) m are the Pochhammer symbols introduced 
earlier. Writing Eq. (|CJ8|) wouldn't normally be too use- 
ful (because Hypergeometric functions are hard to eval- 
uate numerically), except that one can use the identity 21 

F[a, b; c; z] = (1 - z) c " a " b F[c - a,c- b;c;z}. (CIO) 

This identity is extremely helpful because the second 
variable in the Hypergeometric function on the right 
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hand side is a non-positive integer for n = 1, 2, 3, 4, and 
so, because of the definition of the Pochhammer symbol, 
the infinite sum in the definition of the Hypergeometric 
function becomes truncated to n terms. We thus obtain, 
for n = 1,2,3,4, 



i n (y) 



k 2 + y 



' exp(--tan- 1 |)A l (Cll) 



and 
Jn(y) 

where 

D n = 
and 

E n = 



1 2 y 
—2- — j ex P ( taI1 ~ ( ti . 1 ))- E "' 

i) + v y 2k + 1 

(C12) 



(2fc + l) 2 +y 



jk-w^ F[l -? l ~ n]2] ^ (C13) 



n! z 2iy 

(2k + l-iy) n ~ 1 L ~ y ' ~ n ' ' 2fc + l+iy J ' 

With these definitions, straightforward evaluation gives 



D : = 1 



^2 = 

D 3 = 
D 4 = 



2(k - 1) 
k 2 + y 2 

4(Ar — l)(2fc — 1) 2 



(fc 2 + j/ 2 ) 2 k 2 + y 2 

8(k- l)(2fc — l)(3fc- 1) 8(3fe — 2) 



(fc 2 + y 2 ) s 



(fc 2 + y 2 )< 



and 



Si = 1 

E 2 = 



Ik 



(2fc + l) 2 +y 2 



(C16) 



Summarizing, we have 



h 



1 ^ , 2 tan~ l2/ ) 
k 2 + y 2 y k 



2(fc-l) ,2 _ x y. 

l2 = W T7r exp{ 'y t ^ fc } 

(4(fc - l)(2fc - 1) 2 1 ( 2 _ lV , 

(fc 2 +y 2 )3 (fc 2 + y 2 ) 2 / ex Pi~- tan fcJ 

l)(2 fc -l)(3 fc -l) 8(3 fc -2) | 2^, 
-y 2 ) 4 (fc 2 + y 2 ) 3 J v y fc 

^2fc + r ; 

4fc 



(fc 2 



Jl = 



J 2 - 



(2fc + l) ; 



((2fc- 



^ 2 -1/2/ 

^exp( tan ( — — 

> 2 + y 2 v v 2fc + 



?7 



- F _^exp(-^tan- 1 (^ I )) 



(C17) 



(C15) and these are to be substituted into Eq. ([22]) . 
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